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We report on our results from a fully relativistic simulation of the quenched bottomonium spectrum. Using an 
anisotropic formulation of Lattice QCD, we were able to retain a very fine resolution into the temporal direction 
for a range of different lattice spacings. At fixed renormalized anisotropy we study the scaling properties of the 
spectrum and compare our results with non-relativistic calculations. 



1. INTRODUCTION 



2. ANISOTROPIC LATTICES 



The non-perturbative study of heavy quark sys- 
tems is comphcated by the large separation of mo- 
mentum scales which are difficult to accommo- 
date on conventional isotropic lattices. Several 
approximations to relativistic QCD have been 
proposed to describe accurately the low energetic 
phenomenology of heavy quarkonia [QJ^. How- 
ever, the (non-)perturbative control of system- 
atic errors in those approximations is very dif- 
ficult and in practise one has to rely on addi- 
tional approximations. The high precision results 
for the spin structure in non-relativistic bottomo- 
nium calculations have been hampered by large 
systematic errors which are difficult to control 
within this effective theory. Higher order radia- 
tive and relativistic corrections are sizeable for 
bottomonium and even more so for charmonium. 
Still more cumbersome are large observed scaling 
violations j^] that cannot be controlled by taking 
the continuum limit. 

We take this as our motivation to study heavy 
quarkonia on anisotropic lattices in a fully rela- 
tivistic framework. This approach has previously 
been used to calculate the charmonium spectrum 
with unprecedented accuracy ^ . Here we ex- 
tent the method to even heavier quarks and focus 
on our new results for bottomonium. Prelimi- 
nary results were reported in |^ . In Section ^ we 
briefly review our action and simulation parame- 
ters and Section || discusses our results. 
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To study the relativistic propagation of heavy 
quarks it is mandatory to have a fine resolution 
in the temporal lattice direction. To this end we 
employ an anisotropic gluon action: 



(1) 



This is the standard Wilson action written 
in terms of simple plaquettes, Pfii,{x). Here 
(/3,^o) are two bare parameters, which deter- 
mine the gauge coupling and the renormalized 
anisotropy, ^ = as/a*, of the quenched lattice. 
The anisotropic gluon action (|l|) is designed to 
be accurate up to 0{al,a^). 

For the heavy quark propagation in the gluon 
background we used the "anisotropic clover" for- 
mulation as first described in ji,^. The dis- 
cretized form of the continuum Dirac operator, 
Q = iriq + Ip, reads 



Q 



mo + Vs Wiji + vt lUo7o - 
— [cs crakFok + Ct (JkiFa] 
- (a^/2)7^A^ . 



(2) 

The five parameters in Eq. |2| are all related to 
the quark mass, ruq, and the gauge coupling as 
they appear in the continuum action. Their clas- 
sical estimates have been given in ||^. Here we 
chose mo non-perturbatively, such that the rest 
energy of the hadron corresponds to its experi- 
mental value (Af(35f = 9.46 GeV). We also 



fix i^s = 1 and adjust vt non-perturbatively for 
the mesons to obey a relativistic dispersion rela- 
tion ( c(0) = 1 ): 

E\V) = i?2(0)+c2(p) p2 + 0(p4)... . (3) 

For the clover coefficients (c^, ct) we choose 
their classical estimates and augment this pre- 
scription by tadpole improvement. The tadpole 
coefficient has been determined from the average 
link in Landau gauge: = 1/3 (C/^(a:))Landau- 
Any other choice will give the same continuum 
limit, but with this prescription we expect only 
small 0{aa) discretization errors. 

Using the anisotropic giuon action we gen- 
erate several thousand quenched gauge field con- 
figuration for a variety of lattice spacings and 
anisotropics. Every few hundred sweeps we cal- 
culate relativistic quark propagators from the in- 
verse of Eq. ^. This leaves us with several hun- 
dred independent measurements for each set of 
bare parameters. The construction of meson cor- 
relators with definite J^'-' from quark propaga- 
tors is standard Q. Our fundamental bilinears 
are constructed from two 4-spinors, one of 16 F- 
matrices and additional insertions of up to 2 lat- 
tice derivatives: 

M{x) = q{x)V^AJ^^m{x) ■ (4) 

We have further improved those operators to give 
states with different projection onto the ground 
state. Here we follow [Q and employ a com- 
bination of various quark smearings and APE- 
smearing for the gauge links. In practise we find it 
often useful to use gauge-fixed box sources, where 
we can vary the overlap by simply varying the ex- 
tent of the box. This allows us to extract reliably 
both the ground state energies and their excita- 
tions from correlated multi-state fits to several 
channels. We demand consistency for fit results 
from different fit methods and fit ranges. In order 
to call a fit acceptable we require its Q- value to be 
bigger than 0.1. In Table |^ we list the parameters 
of our simulation. 

3. BOTTOMONIUM SPECTRUM 

In Figure we show the continuum limit of 
our results and selected exotic candidates in the 
bottomonium spectrum. 



The later are results from only one lattice spac- 
ing at (6.1,4). To convert our lattice results into 
dimensionful quantities we used the ^Pi — '^Si 
splitting which is ~ 440 MeV experimentally. It 
is well-known that, without dynamical sea quarks 
in the gluon background, the definition of the lat- 
tice spacing is ambiguous and one cannot repro- 
duce all experimental splittings simultaneously. 
For the purpose of this work we admit to the un- 
certainties of the quenched calculation and focus 
on other systematic errors instead. 

As mentioned above we are tuning vt non- 
perturbatively at each mass, so as to reproduce a 
relativistic dispersion relation to within 1 — 2%. 
We tested for systematic errors on spectral quan- 
tities as we vary vt and find that they are small, < 
1%. We have also performed a finite volume study 
and find that volumes with L k, 1 fm are suffi- 
ciently large to accommodate the small bottomo- 
nium ground states. For even smaller volumes we 
could observe some shifts in the spectrum, namely 
the hyperfine splitting is suppressed. 

The most interesting aspect of our work is the 
determination of the spin-structure in bottomo- 
nium from a relativistic formulation. In Fig. ^ 
we show the continuum extrapolation of the hy- 
perfine splitting along with selected results from 
non-relativistic calculations. Clearly there are 
large scaling violations which result in a rather 
large continuum limit with sizeable errors. We 
also notice some deviation from quenched non- 
relativistic calculations at finite lattice spacing. 
We interpret those discrepancies as due to higher 
order relativistic corrections as well as a differ- 

Table 1 

Simulation Parameters. Using ^Pi — ^Si to set 
the scale, we adjust the bare quark mass such that 
the 1 corresponds to the experimental values 



in hb{T). 




N'i X Nt 


configs 


l~- [GeV] 


(5.7, 4) 


8^ X 64 


291 


10.58(67) 


(5.9, 4) 


8^ X 64 


600 


9.25(61) 


(6.1, 4) 


16^ X 96 


660 


10.40(30) 


(6.3, 4) 


16^ X 128 


450 


12.60(67) 


(6.5, 4) 


16^ X 160 


710 


13.04() 


(5.9, 5) 


8^ X 96 


780 


9.57(44) 


(6.1, 5) 


16^ X 128 


370 


10.11(97) 



ent tadpole prescriptions in g. In Refs. 
improved gluon actions were used which can also 
result in different scaling violations. From po- 
tential models and results in we should ex- 
pect dynamical calculations to raise the hyperfine 
splitting even further. We also like to point out 
the very good agreement for our simulations with 
different renormalized anisotropy (4 and 5). This 
indicates that we are working already close to the 
Hamiltonian limit of lattice QCD. 

For the fine structure, Fig. |[ we observe 
a closer agreement with NRQCD calculations. 
However, the large scaling violations result in an 
overestimate of the experimental value by 3a. 

For the first excited states, 25', we also find an 
enlargement compared to experiment and in line 
with previous estimates of this quantity. We feel 
that anisotropic lattices provide a more reliable 
tool to resolve such high excitation, but we leave 
a more detailed analysis to future work. 

In conclusion, we have demonstrated that 
a relativistic simulation of bottom quarks is 
possible on anisotropic lattices. Work is in 
progress to treat also light dynamical quarks in 
an anisotropic fashion. 
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Figure 1. Quenched bottomonium spectrum. We 
plot continuum estimates as full circles. Selected 
results for exotic candidates (hybrids) from finite 
lattice spacing are shown as open squares. The 
^Pi — '^Si splitting is used to set the scale. 
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Figure 2. Hyperfine splitting. We show 
the continuum extrapolation from a linear-plus- 
quadratic ansatz in as. 
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Figure 3. Fine structure. As in Fig. 



